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PERCOLATION 



ELAHE ZOHOORIAN AZAD 

Abstract. In this article, we study a type of a one dimensional percolation 
model whose basic features include a sequential dropping of particles on a 
substrate followed by their transport via a pushing mechanism (see [S. N. 
Majumdar and D. S. Dean, Phys. Rev. Ltt. A 11, 89 (2002)]). Consider 
an empty one dimensional lattice with n empty sites and periodic boundary 
conditions (as a necklace with n rings). Imagine then the particles which 
drop sequentially on this lattice, uniformly at random on one of the n sites. 
Letting a site can settles at most one particle, if a particle drops on an empty 
site, it stick there and otherwise the particle moves according to a symmetric 
random walk until it takes place in the first empty site it meet. We study 
here, the asymptotic behavior of the arrangement of empty sites and of the 
total displacement of all particles as well as the partial displacement of some 
particles. 



1. Introduction 

Fundamental in the domain of percolation is the manipulation of dynamic sets: 
sets that can grow, shrink or otherwise change over time. Some algorithms, like for 
example the Kruskal or Prim algorithms, for the research of the minimal covering 
tree of a graph, involve the grouping of some distinct elements into a collection of 
disjoint sets, and implementing two operations, UNION, that unites two sets, FIND 
that finds which set a given element belongs to, see [CLR90 . 

In a basic model, clusters with different masses change, over time, through space 
and when two clusters are sufficiently close they merge into a single cluster, with a 
probability quantified, in some sense, by a rate kernel R depending on the masses, 
the positions and the velocities of the two clusters. However, such a model, includ- 
ing the spatial distribution of clusters and their velocity, is still too complicated for 
analysis. A first approximation was suggested independently by Marcus |Mar68j 
and Lushnikov |Lus73[ ILus78| , considering kernels depending only on the masses of 
the clusters. 

A Marcus-Lushnikov process )Ald99] is a continuous-time Markov process whose 
state space is the set of partitions of n or, equivalently, the set of measures /i = 
Y^k "t^ ^ k with Sfc kn(k, t) — n, on the set of positive integers N. The fc's stand 
for the sizes of clusters and n(k,t) is the number of clusters with size k at time 
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t. The size-fc clusters provide a fraction fcn ( fc -*) f the total size n. A Marcus- 

1 n 

Lushnikov process evolves by instantaneous jumps according to the rule l each pair 
of clusters (c\,C2) merges at rate R(ci,C2)/n\ which R is the rate kernel of the 
process. 

The sizes of trees in the forest of the spanning-tree model of Yao, or the sizes 
of blocks of cars in the classic parking model |CM04j . form an additive Marcus- 
Lushnikov process in which the rate kernel is R(ci, C2) = c\ + c-i. The sizes of the 
connected components of the random graph of Erdos-Renyi [ER60 , form a mul- 
tiplicative Marcus-Lushnikov process in which the rate kernel is R(c\,C2) — cic 2 . 
The average costs of the Union-Find algorithms, in the model of Erdos-Renyi, were 
studied by Knuth & Schonhage [KS78] . and Stepanov |Ste70| . In two cases, the 
clusters are the connected components of a graph, and the merging of two clusters 
is caused by the addition of an edge between elements of these clusters. We can, in 
this article, suppose that the initial state consists of n cluster of size 1, which corre- 
sponds to a graph completely disconnected with n vertices but any edge. There is 
n — 1 merging between the initial state monodisperse, Si, and the final state, ^S n , 
of the Marcus-Lushnikov process. As it soon will be seen, the model on which we 
work comprises the additive case: the evolution of the sizes of clusters is described 
here by an additive Marcus-Lushnikov process. 

In first analysis, we can distinguish three different regimes in the evolution of 
the additive Marcus-Lushnikov processes, according to the size A% 1 of the largest 
cluster after the A:-th jump, with the interpretations concerning the fragmentation 
of trees |AP981 IPav77j or the analysis of hashing algorithms |CL02| : the sparse 
regime for the case in which if i/n = o(n — k), 1 /n tends to in probability; 
the transition regime, when n — k = 0(y / n), several clusters of size 0(n) coexist, 
and, once renormalized, clusters' sizes converge to the widths of excursions of a 
stochastic processes related to Brownian motion; and finally the almost full regime 
for the case if n — k = o(y / n), x jn tends to 1 in probability, and a unique giant 
cluster of size n — o(n) coexists with smallest clusters with total size o(n). 

2. Main theorems 

Considering a Marcus-Lushnikov process, at the k-th jump, two clusters with 
respective sizes (Sk.n, Sk,n), <Sfc,n > Sk, n are merged, at a cost that may depend on 
the sizes (Sk, n , Sfc,n)- For instance, in some implementations, a label is maintained 
for each element, signaling the set it belongs to, and when merging two sets, one 
has to change the labels of the elements of one of the 2 sets. Yao |Yao76| , Knuth & 
Schonhage |KS78] . studied two algorithms Quick-Find and Quick- Find- Weighted. 
Quick-Find updates the labels of one of the two sets, selected arbitrarily, leading 
to a cumulated cost = X^l'Li Bk,n> m which = Sk, n with probability 1/2 

and Bk : n = Sfe.n with probability 1/2. Quick-Find- Weighted updates the smallest 
set at a cost of Ck, n = s/c,n, leading to a cumulated cost C^^ M/ = X^fcli s k.n- In 
other contexts where coalescence of two sets occurs, costs of interest are Lk, n , the 
size of one of the two sets chosen randomly with a probability that is proportional 
to its size, i.e. L fci „ = S k . n with probability Sk, n / '(Sk,n + s k,n) and L k , n = s fci „ with 
probability s k>n / (S k 
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In our model, where each particle while falling on an occupied site moves ac- 
cording to a symmetric random walk until it finds an empty site, the merging cost 
of two clusters of particles, at the dropping moment of the fc-th particle, is the 
movements of this particle in the cluster on which it falls, until it finds an empty 
site. We indicate this movements by Mk, n > which is also the necessary time for 
fc-th particle to find an empty site. Obviously, Mk. n depends on the size of the 
corresponding cluster (see the definition of the cluster in Section [3]). The partial 
cumulated cost is then 

m 
k=l 

which is interpreted as the total movements of the m first particles (1 < m < n). By 
the two following theorems, we study the concentration of the partial cost C n j an i , 
and the limit law of the total cost C n<n -i, when n tends to infinity. 



Theorem 2.1. For each r\ £ (0, 1), and each e positive, 

C n , \an-\ a 2 (a 2 - 3a + 3) 



liml 



sup 

,ae[0,l-7)] 



6(1 -a) 3 



> e 



Theorem 2.2. We have, 

C n ,n—1 lot 
n 5/2 g 



where £ is a random variable in which the distribution is characterized by its mo- 
ments: 

F(t k \ - n 
l? ' ~ 2(7fc-2)/2 r (5fc-l)" fc ' 

with 

fc-i 

5fe = 2(5fc — 6)(5fc — 4)afc_! + aja^-j k > 2; ax = v^2 . 



The suite of this article is organized in the following way: in Section [3j we 
explain the embedding of the additive Marcus-Lushnikov process in our model, and 
we calculate the probability of the merging of two arbitrary clusters. In Section 
4.1| using only the properties of the symmetric random walk, we can calculate the 



two first moments of the partial cost, M^ n . Theorem 2.1 is proved in Section 4.2 
thanks to the convergence of the additive Marcus-Lushnikov process to the certain 
solutions of the Smoluchowski equation, derived by the analytical arguments in 
|Nor99| : we use, more precisely, Theorem 3.1 of [CM04] . In Section [5] we show that 
the cumulated cost of our model can be approximated by an additive functional on 
Cayley trees, induced by the tolls (?i 2 ) n >i (Proposition [l]), which makes it possible 
to apply the results of [Z 
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FIGURE 1. A sample of tries f(p) and the resulting 3 clusters. 
Here n= 10 = 6 + 2 + 2. 

3. Embedding of the additive Marcus-Lushnikov process 

We start with a description of the additive Marcus-Lushnikov process that helps 

to understand its relations to the parking scheme, generally: at the fc-th step pick a 

I p| 

first cluster P with the probability ^ among the n — k + 1 clusters, and call it the 
'predator' (being a size-biased pick, one obtains, generally, a cluster larger than 
the average cluster); then pick the 'prey' p uniformly among the n — k remaining 
clusters, and suppose that P eat p, producing a unique cluster with size \P\ + 
Remark that if, alternatively both clusters are size-biased picks (resp. if both are 
uniform picks), we obtain the multiplicative Marcus-Lushnikov process (resp. the 
constant kernel Marcus-Lushnikov process, also called Kingman's process) . 

Consider a lattice with n sites at a circle, on which a set V = {1, . . . , n — 1} 
of n — 1 particles drop successively and eventually stick. Letting a site can settles 
only one particle, each particle p drops on a random site f(p). If the first chosen 
site f(p) is on an empty site, the particle stick there. On the other hand, if the site 
f(p) is occupied, the particle executes a symmetric simple random walk and finally 
it stick on the first empty site which it meets. The first chosen sites (f(p)) peV arc 
assumed independent and uniform on the n sites, numbered from 1 to n. 

In this model, the clusters are formed by the occupied sites, with the following 
conventions: 

• there are as many cluster as there are empty sites, 

• a cluster contains an empty site and the set of consecutive occupied sites 
before (going clockwise) this empty site, 

• the size of the cluster is the number of sites constituting it, including the 
empty site, 

• if an empty site follows another empty site, it is considered as a size-1 
cluster of its own. 

The initial configuration, with n empty sites, has thus n size-1 clusters (the 
monodisperse configuration). Each time that a particle sticks, two clusters merge, 
with conservation of the mass, as the empty site that disappears and the particle 
that replaces it both count for one mass-unit. The final configuration, once the 
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n — 1 particles are sticken, is constituted of a unique cluster with size n, and the 
unique empty site, uniformly distributed on {1,2, ... ,n}. It turns out that the 
sizes of clusters form an additive Marcus-Lushnikov process, with kernel K(x, y) = 
(x + y)/n: 

Lemma 1. Give that k particles already sticken (that £ = n — k sites are empty), 
and consider then two clusters with sizes x and y. The probability that these two 
clusters merge at the next drop, p n ,k(x,y), is 



(1) X + V 



n{n — k — 1) 

Proof. Let b\, b 2 , . . . , b h ; 1 < h < k — 1 design the non empty clusters just before 
the k-th particle drops (by convention, b\ design the cluster which contains the first 
dropped particle). As the k—1 first particles choose uniformly their sites, the order 
of these h clusters, from bi on, is a random uniform permutation. On the other 
side, it is not very hard to see that the empty clusters (the clusters with size 1) arc 
merged uniformly on all their configurations: 

Let h, I2, ■ ■ ■ , lh being the number of clusters of size 1 separating respectively the 
non empty clusters (i.e. the clusters of size more than 1) 61, 62, • • ■ , bh- Thus 
dropping the fc-th particle, two empty clusters merge conditioning that these two 
clusters be contiguous and the particle drops on the first cluster (clockwise), in 
other words, conditioning that the fc-th particle drops on one of the 



E ft-^H 

ie{l,...,h} 



empty sites surrounded itself by two empty sites, on the right and on the left. The 
conditional probability that two empty clusters merge, knowing the position of the 
clusters, is thus 



n 

ie{l,...,h} 



This probability does not depend on the sizes of the nonempty clusters, but only 
on their number h, and the position of the r empty clusters among the h nonempty 
clusters (note that k — 1 + h + r = n). Let us pose now r := J2ie{i h} * ne ^ un 
number of empty clusters, and let us note C := {(h,h, ■ ■ , h); Sie{i h} = r l' 
the set of all configurations of ■ ■ ■ , lh), often called compositions of r with h 

pieces. It is well-known that 



CI = Card C = 



r + . 
h 
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Then the conditional probability that two empty clusters merge at the fc-th drop, 
knowing the number of empty and nonempty clusters, is 

YjC Si£{l,...,h}('i ~ !) + 



Pr.hj, 



1 



\ " 



nC hc 



[r-h+ Card {i; h = 0}] 



n(r - 



C r h (r-h)+ Yl 

o<e<h- 



-h+t 



\2J_ 
h 



Remark that J2o<i<h-i ^h-£ +i it) can be interpreted as the number of compo- 
sitions of r in h pieces, such that a null piece be marked, or underlined, whereas 
C^-Ji can be interpreted as the number of compositions of r in h — 1 pieces, such 
that one of the h interstices between the pieces be marked, or underlined: we obtain 
a bijective correspondence between two sets inserting one zero additional into the 
site of the underlined interstice, and underlining the zero so inserted. In addition, 
in the additive Marcus-Luslmikov model, the probability that two clusters of size 1 
merge, at the stage k, if there is r pieces of size 1 and thus h = n — fc + 1 — r pieces 
of sizes higher than 1, is also 

'r\ 2 
s 2 J n(n — k) 
under the terms of Lemma [TJ 

Now let us consider the merging probability of two clusters of respective sizes 
x > 2 and y > 2. Let us note N XjV the number of empty sites met while going 
clockwise, from the cluster of size x to the cluster of size y: N x>y is uniform on 
{1,2,..., n — k}. Obviously, if N x , y ^ {1, n—k}, the two clusters are not contiguous, 
and cannot merge dropping of the fc-th particle. For N x<y € {l,n— k}, let us note 5 
the exit direction of cluster, + or — according to whether the particle leaves there 
in the clockwise direction or in the opposite direction. Let us note r the size of the 
cluster in which the particle drops. We have then, 



Pn,k{x,y) 



rh 



And as (see Section 4.1 1 



r) 



2t 

T - 1 

2t 



if 5 = 



if 5 = 



we obtain well p n ,k{x,y) 



The merging probability of a cluster of size 



x > 2 with one of the r clusters of size 1 in this model, namely 

£+1 

n CI ' 
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coincide also with the probability in the additive Marcus-Lushnikov model, namely 

r(x+l) ' n 

n(n—k) ' 



From now, assume the Marcus-Lushnikov process to be embedded in a drop 
particle scheme. In particular, we preserve the interpretation of Lj >n (resp. Rj, n ) 
as size of the cluster which is chosen by the j-th particle (resp. cluster which 
merges with the preceding cluster when the j-th particle sticks). We indicate by 
Pmle conditional probability that Lj >n (which we will interpret as the size of the j- 
th predator before his meal) is equal to k when the cluster created by the dropping 
of j-th particle (the j-th predator after its meal) is of size to. According to the 
asymptotic behavior of p^ifi) ■, when m is large, we hope to reach a certain intuition 
of the respective values of Lj t „ and Rj, n - It proves, for combinative reasons, that 
Prnk d° n °t depend on j or n: we have, for example, 

Pm,k —P m ,k - r l L m-l,m - «J 

= P(-R m _i, TO = m-k) . 
In what follows, we shall remove thus the exponent of plfk • 



Lemma 2. 



k k - l (m - k) rn - k - 1 (2k - 1) 



k J 4(m — l)m r 



-l 



Proof. Let us calculate the probability that the j-th merge utilizes a cluster 

of size k and a cluster of size m — k, knowing that the result of the merge is of 
size to: as it is a probability concerning the evolution of the sizes of the clusters, 
it is the same one for all the models where the evolution of these sizes is described 
by an additive Marcus-Lushnikov process. It is thus enough to calculate q£fl , as 
that was done in the parking model by Chassaing and Marchand in CM04J. Here, 
we point out this calculation for the convenience of the reader. Among the n J 
configurations for the j first drops, there is 

'i-i' 



2 



n m m -\n - my m (n - j - 1) 



configurations in which the j-th drop form a cluster of size m: there is (^~_} 2 ) 
choice for the m — 2 other particles forming the cluster of size m, n positions for 
this cluster, and once the position and the particles are chosen, there is m m ~ 2 ways 
to build this cluster of size m. The j — m + 1 other particles can be sticken of 
(n — m) J ~ m (n — j — 1) ways on the n — m — 1 sites which are reserved to them. 
Among the configurations in which the j-th drop form a cluster of size m, there 

is 



i-i 

1 , m — k — 1 



i,\m-t-2 



nk (m — k) 



x(k + to — k)(n — my m (n — j — l) 

configurations where the predator is of size k: there is (fe_ 1 J „7-fe-i) ch°i ce for the 
particles of the two clusters intended to be merged, n positions for this set of two 
adjacent clusters, and once the particles of the two clusters and the position are 
chosen, the j — 1 first particles can stick in k k ~ 2 (m — k) m ~ k ~ 2 (n — my~ m (n — j — 1) 
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ways. This calculation holds if we place the cluster of size k initially, and there is 
then k choice for the site where the j-th particle drops. This calculation holds also 
if we place the cluster of size m — k initially, and there is then m — k choice for the 
site where drops j-th particle. It leads to 



k J (m — l)m m 2 

In our parking model, let x denote the probability that there exist two clusters 
of size k and m — k side by side, at the dropping moment of the j-th particle. Then 
the probability that the cluster of size k is before (resp. after) the cluster of size 
m — k is x/2. If the cluster of size k is on the left, the predator is of size k if the 
j-th particle falls on the one of its k — 1 occupied sites and exit from the right, with 
the probability 

x(k - 1) 
4^ ' 

or if the j-th particle falls on the single empty site of the cluster of size k, with the 
probability 

x 

2n" 

If the cluster of size k is on the right, the predator is of size k if the j-th particle 
falls on one of its k — 1 occupied sites and exit from the left, with the probability 

x(k - 1) 
4n 

So the cluster of size k is the predator with the probability 

xk 
2n' 

and the cluster of size m — k is the predator with the probability 

x(m — k) 



We deduced that 



and that 



2n 

xrn _ (j )7l ) 
2n ~ Qm > k ' 

* in , k " m, k 



Finally 



Vra.k 



'm - 1\ k k ~ x {m - k) m - k - 1 
k-l) (m~l)m m - 2 ' 

as expected. □ 



The Lemma [2] and the Stirling formula entails at once that 
Corollary 1. 

(2) Vfc > 1, lim p m ,m~k = 



m— >oo 
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Thus the limiting distribution of the size of the last prey is the Borel distribution, 
in particular related to the explicit solutions of Smoluchowski equations [Ald99] , and 
to the function of tree or Lambert function [Knu98]. Thus, in law, Rm-i.m = O (1). 
However, note that this distribution has an infinite expectation, which is coherent 
with the fact that E [R m -i, m ] = © What we retain of these calculations, is 

that provided L^ >n + Rk, n be large, Rk,n or Sk, n be negligible compared to Lk. n - 



4. Partial costs 

In this section, firstly, we calculate the first and the second moments of the 
movements of the fc-th dropped particle, M^ ni which we will need for the proof of 
Theorem |2.2| Then, we deal with the first part of the demonstration of Theorem 
2.1 which is in fact a corollary of Theorem 3.1 of CM04 , stated here as Theorem 
4.1 



4.1. Moments. As in the previous section, an arbitrary cluster of size s consists of 
s — 1 particles sticken successively and an empty site in the s-th position. Consider 
a particle drops in this cluster, on the one of these s sites chosen randomly (in the 
uniform way). This first choice, noted Xq, is thus a uniform random variable on 
{1,2,..., s}. Consequently, 

»=i 



and 



E(* 2 ) = ^i 2 / S 



,;2 



(* + l)(2* + l) 



Consider now the variable 

h 

X h = Xo + J2 Yi ' ft = 0,1, 2,..., 

i=l 

representing the position of the particle in the cluster of size s, after h step. The 
Yi are the Bernoulli random variables of parameter | with value in {—1,1}. We 
indicate by D s the number of steps inside the cluster of size s, before the particle 
sticks (that one can also see as the time of receive to the edge of the cluster). We 
have then 

¥(X Ds = s) = F(X Ds = s\X = s)¥(X = s) 

+ P(X Ds = s\X f s)F{X ? s) 

1 1 s - 1 

s + 2 s 
■s + l 



2s 

Consequently, F(Xr> s = 0) = ^j". Remark that indicates the last site before 
the cluster (in the clockwise direction), site which is empty. 
The processes Mh = Xfi — h and 

M' h = X% - 2(3ft - 2)Xl + h(3h - 1), 
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are martingales. By stopping theorem, 

E(X* 3 -D a ) = E{X 2 -0), 

which gives 

(3) E(D S ) : S " " ' 



6 

Since the variables Xd 3 and D s are not independents, the calculate of the second 
moment of D s starting from A4' h is not direct. To circumvent the difficulty, we add 
the site to the cluster of size s, thus obtaining a cluster of size s + 1, such that 
we fall in a symmetrical situation: we consider then an initial position Xq uniform 
on{0,l,...,s}. We define then D s and X^ in a way similar to D s and Xp s , but 
these two new variables are now independents. We have thus 

E(D 2 ) = E (D 2 s \X = 0)^- 

s + 1 

+ E(D 2 \X Q ^0)^- 
s + 1 

(4) = ~^-E{Dl). 



Moreover, 



E(X 2 ) 
E(A" 4 ) 



s(2s + l) 
6 ' 
6s 5 + 15s 4 + 10.s 3 



= «) = o = P (*Ds = °) 



30(s + l) 

1 

2 

Stopping theorem for the martingale 

M h = Xh~ 2 ( 3 ^ ~ 2 ) x h + h{3h - 1), 

gives then 

E[X^ - 2(3D S - 2)X 2 bs + D S (3D S - 1)] = E(X* + X 2 ), 
from which it is deduced that 

, 2 ._ (,s 2 -l)(3 S 2 -7) 



(5) E(D 2 S 



45 



4.2. After [an]— th drop. Here, Theorem 4.1 gives the expression, in terms of 
the solution q(k,t) of Smoluchowski equation, of the limit function </s s (a) for the 
partial cost 

[cm] 

^n,\an~\ = ^ y ^ {Sk,n , S k ,n , Uk,n) j 
fe=l 

once C ? r i is normalized by n. This theorem covers a wide class of costs, because 

n, | cm \ J ' 

the general expression q (Sk,m Sk, n , Uf., n ) of the instantaneous cost of the k-th jump 
utilizes a randomization parameter Uf., n > uniform on [0, 1]. The asymptotic behavior 
of the partial cost is expressed according to the conditional instantaneous cost 

<;{x,y) = E [s{S k>n ,s ktn ,U ktn )\ (S kn) s k , n ) = (x,y)} . 
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Theorem |4.1| requires a hypothesis little restrictive of polynomial growth of the 
moment of order 2 of the instantaneous conditional cost, 

(6) Vx,y€N, h(x,y)= f s{x,y,u) 2 du<Ax n y m , 

Jo 

for A, m and n well selected^] The cost <f is supposed nonnegative, and (Uk, n )keti.neN 
denote a sequence of independent random variables uniformly distributed on [0, 1]. 
We note ^ the increasing function of [0, 1) in M + defined by 

<p'(a) = r'^^E^.O (^) ?(MMM) dt, 
•'° ken ten 



and 

q(k,t)= [Kl ~ e ^ ] exp(-fc(l-e-*)) 



We have thus 

Theorem 4.1 ( |CM04p . For all rj > 0, 

ci 



sup 

ae [0,1-7)] 



, [cm] 







We pose now some notations concerning our model: p(k) denotes the particle 
concerned with the fc-th jump, the one which verifies 

# {p 1 1 < P < n - 1 and T v < T p(fe) } = k, 

(where T p indicates the moment when the particle p sticks). Let us note fa(p(k)), 
or fo(k) to be short, the first try of p(k), and note fj(k),j > 1, the j-th try of p(k) 
at the time of its search for an empty site. Let us note Hk la cr-algebra generated by 
the trajectories (fj{t))j >Q i< t <fc_i; the particles p(i); 1 < i < k — 1, and by /o(fc). 
Conditioning by "H^, we obtain, according to the relations ^ and ([5]), 

L 2 - 1 

(7) E[M k>n \H k } = ^"- , 



and 

(8) E [MlJH k ] = (Lln ~ 1} 4 ( 5 3LL " ?) • 

The demonstration of Lemma [2] revealed that L k>n can be written 

L>k,n = Sk ill.. , s k n + St n l S k , 



4.1 



we must 



where Vfe lTi indicates randomly a number in [0, 1]. To apply Theorem 
write Mk t n in the form 

For that we must draw randomly L k , n in the set {sfc,n, <Sfc,n}, using as 
explained above, then we must randomly draw the first test of the fc-th particle 
among the Lk tU sites of the clusters in which it falls, for example in the form 
\Lk jn Wk jn \, where W k , n indicates randomly another number in [0, 1], independent 
of Vk,n- Finally, it is necessary to simulate the random walk of the fc-th particle, 

4n jHj, N denotes the set of strictly positive entire numbers. 
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for example using a sequence (Y kinj e) e>1 of independent random variables {±1} 
symmetrical. It can be done, in a traditional way, by using the coefficients of 
dyadic expansion 

dk,n,l 



to reconstitute 



Uk,n — X! )( 



£>1 

, like below 

Edk,n,2t-1 

e>i 



dk,n, 



4£-2 



2 £ 



We have thus 



£>1 

Yk,n,l — ^dk.nAi — 1 • 

<;{x,y) = E[Mk, n \{Sk,n,Sk,n) = (x,y)] 

= E[E[M ktn \H k ]\(S k!n ,s k , n ) = (x,y)} 

1 

6 
f 

6 

f / x 3 + y 3 
6 \ x + y 



E [ ( L k,n ~ l) | (Sk,n, S k ,n) = (x, y)] 



E 



x 2 l 



1 



In the same way, 
h{x,y) 



E [Ml n \ (S k , n ,s k , n ) = {x,y)] 



E 



(iL„-l)(3i?„-7) (£*,„,**,„) = 



45 

(x 2 -l)(3x 2 -7)x + (y2-l)(3y 2 -7)y 
45(x + y) 

such that ft be a polynomial and satisfies thus the assumption We can then 
apply Theorem 4.1 to ^k,n and <>(x,y) . Thus, for all ij > 0, 



sup 

a£[0,l-7j] 



a 



n, \an~] 



with 



y s {a) = 
1 

12 Jo 

_ 1 
_ 6 



YjYI ( fc3 + ;3 ~ k - ^ ^ dt 

km zeN 

( 9 (.,t), a; 3 - a; ) (<?(., t),l) dt . 



To finish calculation, we use the fact that q(.,t) is expressed in term of the Borcl 
distribution, or of the Lambert function T (cf. |Jan01| ). More precisely, if B a 
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designs a Borel random variable of parameter a, < a < 1, we have, for k > 1 : 



'{B a = k)= { ^^e~ k \ 



and, for a < 1, 

e \B r ; 



1-a 



a-«) 3 ' 

„ r 2a + 1 

E = (1^)5 " 

However it is noticed that, for a = 1 — e _t , 

g(M) = (l-o)P(B„ = fc), 
(q(.,i),:r fe ) - e-*E[S;], 
(<?(.,t),x 3 -a;) = 3e 4t -2e 3t -l, 

from where the calculation of ^(a) carried out higher, which leads to 

c , . a 2 (a 2 -3a + 3) 
^ (a) = 6(1 -a)3 ' 

5. Cumulated cost 

To obtain Theorem |2.2| it remains to estimate the error made approximating 
Mk,n by {sk,n + Sk. n ) 2 (Proposition!!]) , which makes it possible to reveal the relation 
between C n>n -i and the additive functional with penalties (n 2 ) n >o, studied in ZAJ. 
By keeping this goal in memory, we demonstrate initially following lemmas. We 
put 

n-l 

L n = ^ ifc,n ■ 
fc=l 

Lemma 3. ||6C n)n _i — L n \\ 2 — o {n 5 ^ 2 ^j . 

Proof. By developing (6C„ )n -i — L n + n — l) 2 , we obtain: 
||6C n , n _i -L„-n+ 1|| 2 = Si + 5 2 , 



where 



and 



n-l 

fe=l 



~ 2 = 2 £ E [(6M i>n - L\ n + 1) (6Aff jn - L 2 „ + 1)] . 

l<i<j'<n-l 

In consequence of the relation Q, for i < j, 

E [E [(6Mi, n - Ll „ + 1) (6M,- „ - L 2 „ + l) | ]] = 0, 



14 



ELAHE ZOHOORIAN AZAD 



thus S 2 disappears. According to (18]), we have also 



E 



(6M k 



1 



Li 



• ^k, 



6L 



23 

IT 



Therefor, 



fc=i 



< 6n / E 



(i 



[an] , 



da 



According to [Pit87j . for < a < 1, {B? an -, 1 /n) n& ^ {B^ ± denotes the size of the 
longest cluster after the fc-th drop), converges in probability to 0, thus 

4" 



limE 



L 



[an] , 



0, 



and Lebesgue's dominated convergence Theorem completes the proof. 



□ 



Lemma 4. 

Proof. 



^^(ife.n + Rk,n) 2 — L n 



k=l 



= O(n 2 \0gn) . 



y^(£fc,n + Rk.n) 2 — L r 



k=l 



fc=l 
n-1 

+ 2 E [•Rfc.nifc.T,] 

Thanks to the Lemmas 4.2 and 4.8 of |CM04j . for all fc G {1, . . . , n 



1} we have 



(9) 

and 
(10) 

Therefor, 



=G(n 2 logn) 



fc=l 



E [ik,„|L* 



n — k 



n—l ~ 

V E [i? fc . n L fc , n ] = / 



o n — [an] 
l 



E 



' [an] , 

n 



da 



E 



L? 

[an] ,n 



n(n — fo:n~| ) 



da . 



Since for < a < 1, (-B? 1 /n) n eN converges in probability to f |Pit87] L we have 



limE 



and 



limE 



L 



Lr 



[an] ,n 



[an] ,n 



= o, 



?i(rt — fan]) 
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and Lebesgue's dominated convergence Theorem completes the proof. □ 

The two last lemmas involve the following proposition: 
Proposition 1. 



= o I n 5 / 2 



,n + Rk 

k=l 

Now, J2k l =i(Lk,n + Rk.n) 2 is precisely the additive functional on the Cayley trees 
induced by the penalties (n 2 ) n >o studied in ZAJ, to which the reader is referred 
for more details. Though, we represent here Theorem 1.1 of jZAj as the following 
Proposition: 

Proposition 2. Let X n be the additive functional defined on the Cayley trees, 
induced by the toll (n 2 ) n >o- Then, 

n- b ' 2 X n A V2 C, 



where £ is a random variable whose distribution is characterized by its moments. 

2(7fc-2)/2 r (S 

where 

k-l 

at = 2(5fc — 6)(5fc — 4)5fe_i + a,jak~j', k > 2 , a\ = 

3=1 



Finally Theorem 2.2 rises from Propositions [T] and [2] thanks to the following 
theorem [ Bil95| : 

Theorem 5.1. Be (X n ) nl =ft, (Y n ) n ^ and X a random variable, such as for all n, 
X n and Y n are defined on the same probability space. If (X n ) n ^ converges in law 
to X and if (\\X n — Y n \\) ne iq converges in probability to then (Y n ) ne fs converges 
in law to X. 
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